The electroweak vacua, collider phenomenology and possible connection with dark 

energy 



Eric Greenwood, Evan Halstcad, Robert Poltis, and Dejan Stojkovic 
HEPCOS, Department of Physics, SUNY at Buffalo, Buffalo, NY 14260-1500 

Higher dimensional non-renormalizable operators may modify the Standard Model Higgs poten- 
tial in many interesting ways. Here, we consider the appearance of a second vacuum which may 
play an important role in cosmology. For the certain range of parameters, the usual second order 
electroweak phase transition is followed by a first order phase transition that may drive the late 
time accelerated expansion of the universe. Such a potential contains kink-like solutions which in 
turn can play a crucial role in reconstructing the global shape of the potential in colliders, as we 
explicitly demonstrate. 



I. INTRODUCTION 

Today, observational data seems to indicate that our 
universe is going through a period of accelerated expan- 
sion. It remains a mystery, however, what is the driving 
force behind this acceleration. This problem is known as 
the dark energy problem. Observationally, the equation 
of state for the universe is w ks —1, which corresponds 
to a constant, or nearly constant, energy density. The 
minimal solution is the true vacuum energy density or 
cosmological constant. If this is indeed the case, this 
may represent the worst discrepancy between theory and 
observation. The value needed to explain the observed 
acceleration of the universe is (10~ 3 eV) 4 . This value is 
some 124 orders of magnitude smaller than the generic 
predicted value (10 19 GeV) 4 . Since our universe certainly 
does not have the generic value of the vacuum energy den- 
sity, there is a serious question why is it so. This problem 
is known as the cosmological constant problem, and we 
would like to distinguish it from the dark energy problem. 
In particular, it is possible to construct models with the 
scalar field where the driving force behind acceleration is 
the scalar field vacuum energy density. Most of the mod- 
els in the literature, including quintessence,k-essence and 
ghost-condensate, are such models 

0,II&a@. While 
these models can provide the mechanism for accelera- 
tion and resolve some issues connected with it, they can 
not solve the cosmological constant problem. In fact, no 
model dealing solely with the scalar field without address- 
ing gravity can solve the cosmological constant problem. 
In the heart of the cosmological constant problem is the 
fact that the matter field Lagrangian is invariant under 
a shift by a constant. Such a shift does not change the 
equations of motion. However, gravity breaks that sym- 
metry. Therefore, it is very unlikely that the cosmologi- 
cal problem will be solved without addressing gravity in 
a fundamental way 0]. 

Here, we will adopt a common route. We will not 
try to solve the cosmological constant problem. We will 
simply propose a mechanism that can explain accelerated 
expansion of the universe using a scalar field, assuming 
that the cosmological constant problem is solved. The 
scalar field that is driving the acceleration does not have 
to be decoupled from the rest of the universe [1, H, 



[ill . Il2j . It can, in fact, be the Standard Model Higgs 
field. The Higgs field is coupled to the other Standard 
Model particles which means that we can test the model 
in colliders. 

The properties of the Standard Model Higgs potential 
are well known. They depend on the Higgs mass and 
self-couplings. If we limit ourselves only to dimension- 
four renormalizablc operators then the electroweak phase 
transition is second order [35j j . However, it is clear that 
the Standard Model is only an effective low energy the- 
ory and there is no need to include only dimension- 
four operators. Inclusion of the higher dimensional non- 
renormalizable operators may introduce very interest- 
ing features. In particular, adding dimension-six and 
dimension-eight operators can introduce two more sym- 
metric minima. Then, for the certain range of parame- 
ters, the usual second order electroweak phase transition 
is followed by a first order phase transition. This subse- 
quent phase transition may drive the late time acceler- 
ated expansion of the universe and yet leave the imprint 
in colliders. Here, we study such a scenario in detail. 



II. MODEL 

We consider the Standard Model Lagrangian for the 
Higgs field, invariant under the electroweak transforma- 
tions. We include the non-renormalizable dimension-six 
and dimension-eight operators 

V($) = -^ 2 $ t $ + Ai(<J> t $) 2 -A 2 ($ t <I , ) 3 + A3($ t $) 4 + Vo 

(1) 

where $ is the standard electroweak Higgs doublet. Vq 
is a constant which is an overall shift of the potential. 

At zero temperature the CP-even scalar state can be 
expanded in terms of its zero-temperature vacuum ex- 
pectation value v and the physical Higgs boson H : 
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The potential as a function of 4> is given by 
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Substituting Eq. ([2]) into Eq. Q and reading off the co- 
efficient in front of the quadratic terms in H, we can find 
the physical Higgs field mass run as 

m H = + 3Aiw 2 + (7/2)A 3 w 6 - (15/4)A 2 v 4 . (4) 

Defining 
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we can then rewrite Eq. ([3]) in a more convenient way, 
V{<f>) = - £o 0?0l0 2 + ^§(0 2 - 0?) 2 (0 2 ~ <&? + Va (6) 



where 
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Here we note that we have restriction on Ai. From 
Eq. ([5]), requiring that both 0i,02 be real, we have the 
requirement that 
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The role of the parameter £0 is to introduce a con- 
trolled fine tuning of V((f>). If £0 = 0, the potential in 
Eq. ([fi]) has two degenerate minima at = 0i and = 02 . 
If £0 7^ 0, the difference between the energy densities of 
the two vacua is 



SV = £ <; 



(9) 



For the sake of definiteness, we take that 02 > 4>i- 

Vo has been added to the potential to specify the 
false vacuum energy density. Wc require that V(0i) ~ 
(10 _3 eV) 4 . This choice represents a vacuum energy den- 
sity that is sufficient to drive the accelerated expansion 
of the universe. We do not explain the appearance of 
such small number. The solution to the true cosmolog- 
ical constant problem may explain it. For example, an 
interesting numerology 10~ 3 eV« (TeV/Mpi)TeV hints 
toward a gravitational origin of this small number. In 
particular, operators suppressed by powers of Mpi might 
be responsible for it. 

V(0) in Eq. ([6]) is the zero-temperature potential. In 
order to study the sequence of phase transitions we need 
to calculate the finite-temperature effective potential. Fi- 
nite temperature effects are approximated by adding a 
thermal mass to the potential. The potential is then 
written as V(<f>, T) = cT 2 2 /2 + F(0, 0), where c is gener- 
ated by the quadratic terms that acquire a 0-dcpcndent 
mass in the high-T expansion of the one-loop thermal 



potential. Note that there are also terms which are pro- 
portional to T 2 4 , however these terms only lead to small 
corrections to the potential [l3| . 

If the mass of a certain species of particles is greater 
than the temperature of the plasma, the thermal correc- 
tions due to this species decouple exponentially. There- 
fore, strictly speaking, one has to multiply the contri- 
bution from each of the species by the step function 
0(T — m) [28]. While this effect may modify the fine 
details of the phase transition the general qualitative fea- 
tures will remain unchanged. 

The effective potential can then be written as 
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Following the procedure in [l4[ and [l5j|, the constant c 
is given by 
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g and g' are the SU(2)l and U(l)y gauge couplings, and 
yt is the top Yukawa coupling. All temperature depen- 
dence is in e(T). 

Now, we can study that change of the shape of the po- 
tential as the universe cools down (as shown in Fig. [I]). At 
very high temperatures, before the electroweak symme- 
try breaking, the potential has a characteristic "U" shape 
with a single minimum at = 0. The whole potential is 
symmetric and we can consider only > semi-plane. 
The vacuum energy density of the Higgs field before the 
electroweak phase transition is generically of the order of 
the characteristic energy scale of the phase transition, i.e. 
~ lOOGcV. However, notice that the Higgs field is sitting 
there only before the electroweak phase transition where 
the temperature of the universe is high and the universe 
is radiation dominated. The Higgs field has zero expecta- 
tion value and the electroweak symmetry is not broken. 
As T falls, the minimum at = becomes a maximum, 
while simultaneously two new minima appear. The Higgs 
field then rolls down the potential to the first minimum 
at = 0i. There, the Higgs field has a non-zero ex- 
pectation value and the electroweak symmetry is broken. 
This is the standard electroweak phase transition which 
is of the second order. Thus, we do not change the stan- 
dard picture of the early universe. We currently live in 
= 0i vacuum, where the vacuum expatiation value of 
the Higgs field is 0i = 246GeV. 

However, the temperature dependent evolution of the 
Higgs field potential does not stop there. The other min- 
imum at = 02 keeps descending as the temperature in 
the universe drops. Eventually, below some critical tem- 
perate, the minimum at = 02 becomes the true global 
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FIG. 1: The temperature dependence of the <j> Higg s field 
potential. As the temperature decreases the minimum at 
cj> = becomes the maximum and the field starts rolling 
down. Simultaneously two new minima appear. The standard 
electroweak (second order) phase transition is over when the 
Higgs field ends up in the first vacuum at (j) = fa. Then the 
second minimum at 4> = 4>2 descends and becomes the true 
vacuum. This drives the first order phase transition that may 
explain the late time accelerated expansion of the universe. 
The picture is symmetric with respect to the vertical axis. 
The values used in plot are: fa = 0.246TeV, fa = 0.8TeV, 
A 3 = 0.154TeV~ 4 and e = 0.015TeV~ 4 . 



minimum, while the minimum at (j) = fa , that we live in, 
becomes a false minimum. 

In order to justify an effective field theory description 
of the non-renormalizable operators, it is important that 
the coefficients with negative mass dimension in the po- 
tential are considerably less than unity (in units of TeV, 
assuming that the new physics comes at the TeV scale), 
as is the case for the values in the plot in Fig. [1] The 
second requirement is that the scalar field vacuum expec- 
tation values are lower than the TeV scale. In the plot in 
Fig. [TJ the value of the second vacuum expectation value 
is fa = 0.8TeV which is getting close to the new physics 
scale. Note that the values chosen for the plot are just 
for the purpose of illustration, more detailed analysis of 
the possible values of the parameters in the potential will 
be done in Section IVl 

The critical temperature is given by e(T) = 0. From 
Eq. [TT]) we get 



^4>lfa\ 
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At T = T c we have SV = and the heights of the two 
minima are equal. For T < T c , the second minimum at 
fa becomes the true minimum with the energy density 
difference between the vacua given by Eq. 

At high temperatures there is also an overall correction 
to the potential that is proportional to NT 4 , where N is 
the number of relativistic particle species in the plasma. 
This contribution could modify the critical temperature 



if there are different numbers of relativistic species in the 
two adjacent vacua [29|, HO, [3l[. We examine this effect 
in Appendix [Cl 

The existence of the lower minimum than the one we 
currently live in indicates that our universe will eventu- 
ally tunnel into the true vacuum. This phase transition 
may drive the late time accelerated expansion of the uni- 
verse. The tunneling rate depends on the energy differ- 
ence between the vacua which in turn depends on the 
parameter e - 

For completeness, we note that one could set fa > fa, 
with fa identified with the standard electroweak vac- 
uum. In this case, the electroweak phase transition will 
be more complicated than usually thought, since it would 
be composed of two subsequent phase transitions before 
the Higgs field settles down into today's vacuum fa. This 
interesting possibility is outside of the scope of this pa- 
per. 



III. TUNNELING RATE 

An important question is how long will the universe 
exist in the false vacuum state 4> = fa? The transition 
from the false to the true vacuum occurs from nucleation 
of bubbles of true vacuum inside the false vacuum. The 
transition probability per unit space-time volume, using 
the semi-classical approximation, is given by 



= Aer SE 



(14) 



where Se is the Euclidean action of the bounce solution 
and A is a dimensionful constant, which depends on the 
loop corrections to the potential Eq. ([3]). However, here 
we are only interested in order of magnitude transition 
rate, thus we will ignore these corrections. To calculate 
Se, we follow the method developed by Coleman fl6j | . 

The one-dimensional Euclidean action per unit volume 
for the tunneling is, to zeroth order in eq 



Si = / dfa VW{fa) 
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where we are again assuming that fa > fa. In the zero- 
temperature limit and thin wall approximation, the ra- 
dius the critical bubble is then 

R = 35i = 3 [^ {fa-faf^l+ifafa+cpD 
° 5V 5 V 8 e (fa + fa)<j>lfa) ■ [ ' 

For an 0(4) symmetric bubble, the Euclidean action is 
then given by 
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At zero temperature, Eq. (14)) gives the decay rate per 
unit volume per unit time. In order for our observable 
universe, whose four- volume is of the order of tjj ubble , 
to remain in the false vacuum, one must require that 



Hubblt 



< 1. Taking t Hu , 
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10 years, we find that 
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sufficient stability for the false vacuum is obtained for 
Se > 400. With the generic value Ai ~ 1, we see that 
vacuum stability requires only Eq < 0.012TeV -4 . This is 
also enough to make the thin wall approximation valid. 

To ensure that the above analysis remains true for the 
early universe, we calculate the temperature dependent 
decay rate in the high temperature limit. At finite tem- 
peratures the (9(4) symmetric bounce is approximately 
the periodic in time 0(3) symmetric solution. In this so- 
lution, the period is 1/T (see 17[). The decay exponent, 
i.e. the Euclidean action, now has the form 



Se = 



T 
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where S3 is the three-dimensional action of an 0(3) bub- 
ble. The new radius of the critical bubble is now 
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The three-dimensional action of the 0(3) bounce solution 
is then given by 



Sa(T) = --ttR(T) 3 5V{T) + 4ttR(T) 2 S 1 
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Therefore the temperature dependent decay rate is given 
by Eq. (HI), Eq. $Z0§ and Eq. (H) as 



exp 
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Te{T)' 



(22) 



where e(T) is given in Eq. (TTJ). We can see here some im- 
portant properties. First, at high temperature, T » T c , 
when (j> = (f>i is the lowest energy state, the transition 
rate is large. Thus most of the universe ends up in that 
minimum. At T = T c , we have that s(T) = 0, thus the 
transitions between the vacua are suppressed. The high 
temperature approximation is valid, at least formally, at 
the temperatures slightly below T c . It is there that the 
decay rate in Eq. ((22)) is maximal and we need to cor- 
rect the zero temperature estimate for eq. Fortunately, a 
slight correction Eq ~ 0.005TcV~ 4 makes the decay rate 
safely small. As shown earlier, for T <C T c , <f> — <j>i is a 



false minimum, but the transition rate to the true vac- 
uum at 4> = 4>2 is suppressed by the bare value of eq. 
We saw that £0 ~ O.OlTcV -4 makes the transition time 
larger than the current Hubble time. In order to incor- 
porate a somewhat stronger constraint for high temper- 
atures, we require Eo ~ 0.005TcV -4 . 

We have to make sure that energetic processes in our 
universe (e.g. cosmic ray collisions) were not able to initi- 
ate the formation of a true vacuum bubble, which would 
in turn encompass most of the visible universe by now. 
Fortunately, it is not a simple thing to create a vacuum 
bubble in a high energy collision. This requires not just 
sufficient energy but a coherent superposition of a large 
number of high energy quanta over a volume large com- 
pared to the characteristic energy. Such processes require 
high densities and high temperatures, not only high en- 
ergies. The height of the barrier between the false and 
true vacua is 
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This is approximately (O.lTeV) 4 for the values <\>\ ~ 
0.246TeV, <\> 2 ~ 0.8TeV and A 3 - 0.154TeV^ 4 . Such 
temperatures are unlikely to soon be achieved in colliders, 
and are probably not achieved over large enough volume 
even in the highest energy cosmic ray collisions. 

The critical bubble radius, Eq. (Ill), is lOTTeV" 1 for 
our parameters (e - 0.098TeV~ 4 , <§> x - 0.246TeV, <f) 2 ~ 
0.8TeV and A 3 - 0.154TeV- 4 ). Given Eq. (J23J) , this 
suggests that we need approximately 



N, 



quanta 



I k*R\ 



\ 3 



v ma ^j v-j/* > io 9 (i/ £o ) 3 . (24) 



More- 



individual excitations coherently superimposed, 
over, N quanta grows very fast, as Eq 3 , making it extremely 
difficult to create a critical bubble in a high energy colli- 
sion. 

In the case of tunneling from dS to AdS vacuum, there 
is an additional suppression due the different asymptotics 
of these space-times. The size of the critical bubble in 
AdS ends up being larger than one might gcnerically ex- 
pect [HI], which further suppresses the transition. 



IV. FUTURE OF THE UNIVERSE 

In this context we can address the question of the fu- 
ture of our universe. Phase transitions give quite different 
predictions for the future of our universe than the true 
cosmological constant. In the case of the true cosmologi- 
cal constant, accelerated expansion never stops. Acceler- 
ation will slowly drive most of today's visible universe out 
of the cosmological horizon. In a distant enough future, 
the whole visible universe will be a gravitationally bound 
system consisting of only Milky Way and Andromeda 
galaxies [H[20[. 
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However, in the context of phase transitions we have 
several different possibilities. First, if the difference be- 
tween the two vacua is very small (i.e the parameter 
£o ~ 0), the phase transition will never be completed. 
The bubble nucleation rate will be very slow and they 
will never percolate, since the background is expanding 
with acceleration. This scenario is similar to the true 
cosmological constant, except for the possibility to have 
a few bubbles here and there locally. 

The more interesting case is when the parameter e 
takes much less fine tuned values, for example those 
which allow for the phase transition to be completed. 
Phase transitions are violent events and many things will 
be different in the new vacuum after the completion of 
the phase transition. The true vacuum at = 02 clearly 
has a different Higgs field vacuum expectation value that 
the vacuum we currently live in. This means that most of 
the Standard Model particles will have different masses. 
In such a universe it is very difficult to imagine life similar 
to ours due to the well known anthropic reasons. 

Finally, for the perhaps most generic value of £o ~ 
0.01TeV~ 4 , the phase transition is just about to happen, 
in cosmological terms. Since the characteristic scale in 
the Higgs potential is of the order of lOOGeV, the re- 
quirement that the false vacuum at <j> = (f>i is shifted 
up from zero by a tiny amount of I0 _3 eV directly im- 
plies that the true vacuum at = 02 is deeply AdS, i.e. 
has a negative vacuum energy density [2l|, |22|, [23|, [24| . 
The transition from the false to the true vacuum will be 
described by the Coleman-De Luccia instanton [25| . Ac- 
cording to [26[ , any initial instabilities in the AdS space 
will quickly grow and cause the collapse of the whole uni- 
verse into a black hole. In such a scenario it is difficult 
to imagine any life at all. 



V. CONSTRAINING THE VALUES OF 
PARAMETERS IN THE POTENTIAL 

The original potential has four parameters: /u, Ai, A2 
and A3 . It would be interesting to find the possible values 
that do not lead to dangerous exotic vacua. The main 
constraint comes from Eq. I|17p . where we require that 
Se > 400. The requirement that today we live in the 
standard electroweak vacuum is that 0i = 246GeV. From 
Eq. (0,01 can be expressed as <f>\ = 0i (Ai , A2 , A3 ) , which 
puts one constraint on the possible values of parameters. 
From this constraint we can express A3 in terms of Ai 
and A2 . In order to make a useful plot of Se in terms of 
the two parameters, we need to fix one more parameter. 
We do that for /1 by setting the value of the Higgs mass 
(Afniggs = v2/z). Since the exact value of the Higgs mass 
is not known, we will make two plots, one assuming that 
the Higgs mass takes its lowest (experimentally) allowed 
value of 114GeV (Fig. \2§ and the other one assuming 
that the Higgs mass takes the value of 200GcV (Fig. [4|. 

If we set the value of the Higgs mass to 114GcV, from 
Fig. [5] we see that possible ranges of the allowed val- 



ues for the parameter Ai and A2 are < Ai < 0.2, and 
< A2 < lTeV~ 2 . Incorporating these limits into the 
constraint equation 0i = 0i(Ai, A2, A3), we get that a 
possible range of the allowed values for the parameter A3 
is given by < A3 < 10TeV~ 4 , as shown in Fig. [3] 




FIG. 2: A scan over the allowed range of parameters Ai and 
A2 in the original potential (|3]) which do not lead to the phe- 
nomenologically excluded new vacua. The constraint comes 
from the requirement that the Euclidean action Se > 400. 
For this plot, we set the value of the Higgs mass to 114GeV. 




FIG. 3: The allowed range of the parameter A3 in the original 
potential ((3]) which does not lead to the phenomenologically 
excluded new vacua. We use Eq. ([5]l and the requirement that 
4>i — 246GeV to express A3 in terms of Ai and A2. For this 
plot, we set the value of the Higgs mass to 114GeV. 

If we set the value of the Higgs mass to 200GeV, from 
Fig. 0] we see that possible ranges of the allowed val- 
ues for the parameter Ai and A2 are 1 < Ai < 1.5, and 
< A2 < ITcV" 2 . Incorporating these limits into the 
constraint equation 0i = 0i(Ai, A2, A3), we get that a 
possible range of the allowed values for the parameter A3 
is given by < A 3 < 10TcV- 4 , as shown in Fig. [3 As 
mentioned before, the effective field theory description 
of the non-renormalizable operators is valid only for the 
values of A2, A3 < 1 (in units of TcV, assuming that the 
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FIG. 4: A scan over the allowed range of parameters Ai and 
A2 in the original potential Q which do not lead to the phe- 
nomenologically excluded new vacua. The constraint comes 
from the requirement that the Euclidean action Se > 400. 
For this plot, we set the value of the Higgs mass to 200GeV. 



1. 




FIG. 5: The allowed range of the parameter A3 in the original 
potential © which does not lead to the phenomenologically 
excluded new vacua. We use Eq. ([5]) and the requirement that 
(pi = 246GeV to express A3 in terms of Ai and A2. For this 
plot, we set the value of the Higgs mass to 200GeV. 

new physics comes at the TeV scale). 

Having the allowed range of parameters in the poten- 
tial we can estimate the range of values that the second 
vacuum fa can take. From ([5]) we can see that fa can 
take values from 300GcV to well above a TeV. Again, 
the effective field theory description is valid only for the 
values of fa smaller than TeV. 



VI. RECONSTRUCTING THE POTENTIAL IN 
COLLIDERS 

The main feature that distinguishes this type of models 
from other scalar field models (say quintessence) is the 
explicit connection with particle physics and the possi- 
bility of reconstructing the potential in colliders. With 



the inception of the LHC, we hope to be able to achieve 
this there. 

In this model, we explicitly used the (j> 8 potential. How 
do we probe the global structure of the potential? If we 
excite the field locally, only around one vacuum, then we 
can only probe the local structure of the potential. In 
order to probe the global shape of the potential we need 
to excite a solution that extrapolates between the vacua. 
These are non-perturbative kink- like solutions. In [l8l ]. 
the author constructed a recipe for the reconstruction of 
the potential when kink-like solutions are present in the 
theory. The base of the recipe is the inverse scattering 
method. In the direct scattering methods we start from 
the known potential and calculate the eigenfrequencies, 
i.e. the energies of the scattered particles. In the in- 
verse scattering method, we start from the known eigen- 
frequencies (acquired presumably in the scattering exper- 
iments) and calculate the shape of the potential which is 
the cause of the scattering. 

In our context, kink- like solutions which connect the 
two vacua correspond to the bubbles of the true vacuum. 
Inside the bubble we have the true vacuum while outside 
is the false vacuum. The true vacuum is energetically 
favored so the volume term contributes to the pressure 
directed outward. However, the surface tension of the 
bubble tends to contract the bubble. The critical bubble 
is the one which is large enough so that these two forces 
are balanced. We showed above that the production of 
the critical bubble of the true vacuum which is capable 
of expanding is very suppressed. However, production 
of a subcritical bubble which will collapse under its own 
tension should not be severely suppressed. By studying 
production and decay of these subcritical bubbles we can 
learn a lot about the global shape of the potential. 

Bubbles of true vacuum are soliton-like solutions, and 
in zeroth order the production cross section should be 
just the geometrical cross section, i.e. t-Rbubbie where 
-Rtmbbie is the geometrical radius of the bubble. The situ- 
ation is somewhat similar to mini-black hole production 
in high energy collisions. Black holes are gravitational 
solitons and their production cross section was shown to 
be just the geometrical cross section, i.e. nR% H UK]. 
However, there is one significant difference. While in the 
of case mini black hole production, according to the hoop 
conjecture, there are no phase space suppression factors, 
in the case of the bubble the suppression factor must be 
present. A bubble is a coherent superposition of a certain 
number of the scalar field quanta, say N. Therefore the 
suppression factor would likely go as e _entropy which is 
roughly e^ 1 "^ 1 ) ~ e~ N . For a subcritical bubble where 
N is a few, the production may be possible [34 ] ■ For 
a critical bubble with large N, the production is highly 
suppressed. Note also that the threshold for the bubble 
production would be Nm, where m is the mass of the 
scalar field quanta. 

If the bubble of the true vacuum is produced, it is 
unlikely that it will be produced in its ground state. In- 
stead, we expect that it will be produced in a highly 
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excited state. Then the decay of such a bubble will give 
off the eigenvalues of the potential. Since bubbles are 
coherent states of a certain number of the scalar field 
quanta, they would dominantly decay into those scalar 
field quanta (in this case the Higgs field). The eigenval- 
ues of the potential would come from the (reconstructed) 
energy distribution of emitted scalar field quanta. Due 
to the spherically symmetric configuration of the bub- 
ble, one of the main signatures of the bubble production 
would be a spherically symmetric distribution of emitted 
scalar field quanta. For the light Higgs, the main decay 
channel will be through b-quarks, which further decay to 
a c-quark, a lepton (which may serve as a trigger) and 
neutrino. In that case, b-quarks will be produced co- 
piously. Because of the high multiplicity, b-quarks may 
not be very energetic. The main question is then whether 
this signature can be distinguished from the QCD back- 
ground. Two things arc important in this context: first, 
high degree of spherical symmetry, and second, many jet 
events. For example, for five (and likely more) jet events 
the trigger threshold may be lowered to about 50GcV, in 
which case the QCD background can be kept under con- 
trol. The heavier Higgs is much easier to analyze. For 
the heavy Higgs, decay channels including W ± and Z 
gauge bosons becomes significant. Once the ZZ branch- 
ing ratio become significant, the experimental signature 
is much more distinct. The other question is the possi- 
bility of the reconstruction of the original energy of the 
bubble through the decay products. Situations contain- 
ing neutrinos in the final stage are specially inconvenient. 
However, based on earlier experience, say with the top 
quark, we know that something like that is possible. An- 
alyzing other decay products, one may identify the exact 
channel of decay and thus estimate the energy taken by 
neutrinos. 

If a scalar field theory is written in the standard form 



Therefore 



L 



1 



d^d^ - v(<f>) 



(25) 



then the equations of motion can be written in the 
Schrodinger-like form 



dx 2 



dx 2 



1pn{x) = Col^nix). (26) 



For simplicity, the field 4>(x) is a function of only one 
coordinate. <f>o(x) is the (unknown) profile function of 
the kink solution. The task is to determine the poten- 
tial V{4>) knowing the eigenvalues ui n . Though the an- 
swer is not unique, additional theoretical input (e.g. to- 
tal energy of the bubble and some perturbative interac- 
tions) can possibly reduce degeneracies. The zero mode 
ipo = d(j)o /dx which corresponds to the eigenvalue ujo = 
plays a special role. In theories where kink solutions are 
present the Bogomolnyi equation directly relates the zero 
mode to the shape of the potential 



d(f>p{x) 
dx 



(28) 



where x(tf>o) is obtained by inverting <j>o(x) which is in 
turn obtained by integrating the relation ipQ = d<po/dx. 
Thus, by finding the zero mode solution to the Eq. (|2l)|) 
we can reconstruct the shape of the potential. The prob- 
lem is that, in most of the cases, the zero mode is coupled 
to all of the higher modes, and we need to solve a coupled 
set of differential equations. However, as shown in [18l |. 
two of the most important theories containing kink solu- 
tions, i.e. sine-Gordon and 4 theory, have respectively 
one and two eigcnfrequcncics. In these cases we can solve 
the equations analytically. In more complicated cases, it 
is possible that the equations need to be solved numeri- 
cally. 

We now come back to our <fi 8 potential. Since the orig- 
inal potential has four parameters (/x, Ai, A2, A3), we can 
infer that in this case there are four eigenfrequencies: 
Ki, K2, K3, and K4. Following the procedure in [l8j|, in 
the special with a high degree of symmetry between the 
eigenvalues, one can reconstruct the potential (for details 
see Appendix A) 



(29) 



where a is a normalization constant. Here we recognize 
the </> 8 potential. In order to get a desired potential with 
two pairs of non-equivalent minima we need to relax the 
conditions on the eigenvalues. In this more generic case, 
the potential can not be obtained in an analytic form 
since the differential equations must be solved numer- 
ically. The result of the calculations presented in the 
Appendix A is shown in Fig. [6] By comparing the po- 

Reconstructed Potential 




<t> 



FIG. 6: Reconstruction of the cj> S potential using the inverse 
scattering method. The specific values for this plot are v = 1, 
A = —4, K2 — k| =5 and k| — k\ = k| = 3. The plot is 
symmetric with respect to vertical axis. 



±V2^(0 O )- 



(27) 

tcntial in Eq. ^ with the reconstructed one in Fig. [5J 
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we can infer the values of the parameters in our origi- 
nal potential. In particular, we can obtain the value of 
the fine tuning parameter Eq. If the solution to the true 
cosmological constant implies that the true vacuum en- 
ergy density vanishes, then the parameter Eq must give 
for the energy density difference between the vacua in 
Eq. © the value of SV w (1CT 3 cV) 4 . 



VII. CONCLUSION 

The electroweak Standard Model is believed to be 
only a low energy effective theory. For this reason we 
may be allowed to introduce higher dimensional non- 
renormalizablc operators. We expect new physics to kick 
in close to a TeV scale. Here we studied how these higher 
dimensional non-renormalizable operators may modify 
the Standard Model Higgs potential in a manner inter- 
esting for cosmology. We considered the appearance of 
a second vacuum in the Higgs field potential. We cal- 
culated finite temperature corrections and showed that 
the usual second order electroweak phase transition is 
followed by a first order phase transition that may drive 
the late time accelerated expansion of the universe. Such 
a potential contains kink-like solutions which in turn can 
play a crucial role in reconstructing the global shape of 
the potential in colliders. We explicitly demonstrated it 
using the inverse scattering method adopted to studies 
of theories where kink- like solutions are present. 

We addressed the future of our universe in this context 
which is quite different from the future dictated by the 
true cosmological constant. 

Since our model does not address gravity in a funda- 
mental way, it does not solve the cosmological constant 
problem. It simply addresses the dark energy problem as 
all of the other scalar field models do — postulates the 
potential consistent with a given equation of state, which 
in this case is if = — 1. The advantage here is that we do 
not need to postulate the existence of a new scalar field 
completely decoupled from the rest of the universe. As 
an additional bonus in this model, the potential can be 
in principle reconstructed in colliders. 

While we do not explain the appearance of a 
small number 10 -3 an interesting numerology 10 _3 eV« 
(TeV/Mp;)TeV hints toward a gravitational origin of this 
small number. Operators suppressed by powers of Mp; 
might be responsible for it. 
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APPENDIX A: RECONSTRUCTION OF 
POTENTIAL 

We illustrate here the inverse scattering problem, 
i.e. reconstruction of the potential from its eigenvalues. 
Knowing that the 4> 4 potential has two eigenvalues, we 
expect that he (f> s potential has four. We take the four 
eigenvalues of the bound-states to be k±, k%, K3, and K4. 
We will first consider the special case where the solu- 
tion can be found analytically, then we will do numerical 
analysis of a more generic case. 

The translational mode is always the lowest, so we can 
write K4 = 0, since we are using the notation m > Kj+i. 
Following [l8| , we find the potential containing n of the 
highest bound states: 

U n {x) = fl + f n + n 2 n (Al) 

where the function f n (x) satisfies 

/;-/n + ^»-l=«n- ( A2 ) 

Defining f n (x) = —w' n /w n we can rewrite Eq. (|A2[) as 

- + U n -lW n = K^Wn. (A3) 

This equation will have two linearly independent solu- 
tions; however, if we require that U n is even under par- 
ity transformations then we must also require i«„(-i) = 
w n {x). This condition then eliminates one of the linearly 
independent solutions. 

Therefore to find U\ we consider 

Ui = fi +f[ + nl (A4) 

and solve 

- <w" + UqW! = nfwi. (A5) 

Now if we take Uq = k 2 , and define v 2 = Kq — k 2 we can 
then write Eq. (|A5[) as 

- w" + v 2 w x = 0. (A6) 

The solution to this is w± = cosh(j/a;), therefore 

fi(x) = — vtB,vh{vx). (A7) 

Substituting into Eq. (|A4[) yields 

Ui = v 2 [1 - 2sech 2 (ra)] + k\. (A8) 

To find U2 we consider 

U 2 = fl + & + (A9) 

and solve 

- w'2 + Uiw 2 = k 2 2 w 2 . (A10) 

To solve Eq. (|A10[) we use Eq. (|A8[) . Now defining z = vx 
and ft 2 = k\ — k 2 we can then write Eq. (|A10[) as 

^ + (A + 2scch 2 (z))«; 2 = (All) 
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where 



A 



P 2 
v z 



(A12) 



In general the solution for W2 involves the hypcrgeometric 
function, however it is instructive to consider a special 
case. If we take 1 jv 2 = 3, the solution for W2 is W2 = 
cosh 2 (z/:r), therefore 



f2{%) — — 2vt&nh(vx). 
Substituting into Eq. (|A9|) yields 

U 2 = 2v 2 [2 - 3sech 2 (ra)] + k 
To find U3 we consider 



(A13) 
(A14) 



U 3 = fl + & + 4 (A15) 



and solve 



W 3 ' + U2W3 = K3W3. 



(A16) 



To solve Eq. (|A16[) we use Eq. (|A14|) . Now defining z = 
vx and define j 2 = k 2 — k 2 we can then write Eq. (|A16|) 



d 2 w 3 
dz 2 



+ (a + 6sech 2 (z)) w 3 = 



(A17) 



where a = — (4 + j 2 /v 2 ). In general the solution for u>2 
involves the hypergeometric function, however it is more 
instructive to consider a special case. If we take -f 2 ' jv 2 = 
5, the solution for W3 is W3 = cosh 3 (^x), therefore 



fd,{ x ) = — 3^tanh(z/x). 
Substituting into Eq. (j A15|) yields 

U 3 = 3is 2 [3-4scch 2 (i/ir)] 4 
Finally to find U4 we consider 



and solve 



w 4 + U3W4 = K 2 W4. 



(A18) 



(A19) 



(A20) 



(A21) 



To solve Eq. (|A21[) we use Eq. (|A19[) . Now using z = vx 
we can write Eq. (|A21[) as 



d 2 W4 

dz 2 



(p+ 12sech 2 (z)) Wi = 



(A22) 



where p = — (9 + k\/v 2 ). In general the solution for 
u>4 involves the hypergeometric function, however it is 
more instructive to consider a special case. If we take 
k\/v 2 = 7, the solution for W4 is W4 = cosh 4 (ra). 
The profile function is defined as 



t>a = a 



dx 
w N 



(A23) 



where N is the highest eigenvalue, in this case 4, and a is 
a normalization constant. Therefore here we can write, 



— tanh(ra) (sech 2 (vx) + 2) 



(A24) 



For significantly high values of vx this can be approxi- 
mated as 



2 " W X 

— tam\(vx) 



(A25) 



The symmetry breaking potential is then defined from 
the Bogomolnyi equation 



where 



ax 



(A26) 



(A27) 



Therefore we can write the symmetry breaking potential 

as 



V(<h) 



-sech 8 (^x) 



1 



& 4 



4a 2 



(A28) 



where we used Eq. ()A25p . As expected, the potential in 
question is </> 8 . However, the conditions that we imposed 
on the eigenvalues introduced a high level of symmetry 
in the potential. In order to get a desired potential with 
two pairs of non-equivalent minima we need to relax the 
conditions on the eigenvalues. 

Now, we consider a more generic case of the potential 
that we need to reconstruct. In this case we search for 
solutions to our given potential, Eq. (p}. Such as in the 
special case, the potential Ui is given by Eq. (|A4[) . Using 
this potential, we can again solve the differential equation 
for wi(x). The solution for w\(x) is given by wi(x) = 
cosh(ra), which then yields /1, Eq. ([AT]) . From Eq. (|A"8|) 
and Eq. (|A10[> . we can then solve the differential equation 
which has a solution 



w 2 {x) = P 1 v_A (tanh(ra)) + Q{ (tanh(i^)) (A29) 

where A is the same as in the special case above, 
Eq. (|A12[) , and P and Q are the Legendre polynomials of 
the first and second kind, respectively From the parity 
condition, we can see that for the Legendre polynomials 
of the second kind \J — A must be an even integer. For the 
Legendre polynomial of the first kind, the only non-zero 
terms are for \J — A = 0, 1. For \J — A = 0, this does not 
satisfy the parity condition, hence it is not an allowed 
solution. For J — A = 1 . we then have the solution 



1^2(2;) = — yl — tanh(ra) 2 = ±sech(vx). 
This then gives that the function f2 is given by 
/2 = ±v tanh(i/x). 
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This is just T/i; hence when we solve for 103(2;) we will 
just obtain Eq. (|A29|) . We will end up going in circles fol- 
lowing this value; it is therefore more economic to explore 
the Legendre polynomial of the second kind. 

From Eq. (|X29|) we can find the function f 2 to be 



h 



(2 - V=X)i/£# 
— 2^tanh(z/cc). 



^ (taiah{yx))Qi 



A (tanh(i/x)) 

(A30) 



From Eq. (|A30|) . we can then find the potential U 2 using 
Eq. ([AT]) . We use Eq. (|A"30|l to find U 2 , then use this 
to find ui3(x). However, there is no closed form solution 
for W3(x) and we need to use numerical methods from 
here out to find the solution for both 103(2;) and w±{x). 
Using Wi(x), we can then find the reconstructed potential 
V(4>o) which we show in Fig. [6] The specific values that 
give the characteristic shape in question are v = 1, A = 
— 4, «2 — k| = 5 and «§ — k\ = «§ = 3. For convenience, 
we reconstructed the potential in two intervals, the first 
one between <fi — and <fi = 0.8TcV, and the second one 
between </> = 0.8TeV and 4> = 00. We then joined the 
plots into a single plot in Fig. [6] 



APPENDIX B: KINK SOLUTION OF THE FIELD 
<j> BETWEEN THE VACUA fa AND fa 

The potential given in Eq. (|6|) contains a kink-like solu- 
tion interpolating between the vacua fa and fa. Setting 
£0 = in the Bogomolnyi equation we find 



(Bl) 



We evaluate this integral using the method of partial frac- 
tions, 



3l<?2 



^\)x = In 



(B2) 

This gives us x{fa). We then numerically invert it to get 
fax). The result is shown in Fig. [7] 



APPENDIX C: CRITICAL TEMPERATURE 

Several effects may change the numerical value of the 
critical temperature of the phase transition. We first con- 
sider the possibility that the two minima do not contain 
the same number of degrees of freedom. Let us denote 
the number of degrees of freedom in the two vacua fa 
and fa by N% and N 2 respectively. The condition for the 
critical temperature is that the values of the effective po- 
tential pTJ|) evaluated at the two vacua fa and fa are the 
same. From (|10[) , after we add the terms that depend on 
N, we have 




FIG. 7: A kink-like solution interpolating between the vacua 
fa and fa for eo = 0. 



V eff (fa,T) = -[ £o -^J fafa 2 - — AqT 4 + const. 

(CI) 



V e ff{fa,T) = -{s Q -^J 4>\4>l - ^N 2 T* + const. 

(C2) 

where Ni, i = 1,2, represent the contribution over the 
relativistic bosonic, Nb, and fermionic, Nf, spin states, 
i.e. 



N = N b + -N F . 



(C3) 



At T = T c , we have V e fj{fa,T c ) = V ef f(fa,T c ) which 
gives a quartic equation for T c 



(Ni - N 2 )T± + c(4 - fa\)Tl - eofa^^l 
We can now solve Eq. (|C4|) for T c to get 



= (C4) 



^\{fa 2 



(C5) 

We can perform a simple estimate of AiV = Ni — N 2 
if we assume that in the vacuum fa all the particles are 
relativistic, while in the vacuum fa the Higgs boson, W ± 
and Z bosons and top quark masses become larger than 
the temperature of the universe and their contribution 
should be excluded. Therefore AN ~ 15. 

Another consequence of the fact that particles not 
much lighter than the temperature should be decoupled 
from the plasma is that the constant c defined in Eq. (fl2|) 
gets modified. If the temperature during the clcctroweak 
phase transition is of the order of lOOGeV, then it is rea- 
sonable to exclude the Higgs boson, W ± and Z bosons 
and top quark contributions in Eq. (|12[) . This is espe- 
cially true in the second vacuum fa. We first need to 
express the constant c in terms of physical quantities. 
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For this purpose, we eliminate Ai using Eq. ((4]), he. 

Ai = + l^/v 2 + \\ 2 J - 7 -X 2 v\ (C6) 

The constant c now becomes 

+ -k x *" 2 --k^)- < C7 > 

We now exclude the Higgs boson, and Z bosons and 



top quark contributions to get 




where yb is the b-quark Yukawa coupling, which we 
needed to include once the top quark is excluded. In 
turn, the critical temperature defined in Eq. (| 13[) will 
change. While the numerical value of the critical tem- 
perature will change, the qualitative features will remain 
unchanged. 
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